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Abstract 

We have used 106 pb -1 of data collected in pp collisions at y/s = 1.8 TeV by the 
Collider Detector at Fermilab to measure jet angular distributions in events with two 
jets in the final state. The angular distributions agree with next to leading order 
(NLO) predictions of Quantum Chromodynamics (QCD) in all dijet invariant mass 
regions. The data exclude at 95% confidence level (CL) a model of quark substructure 
in which only up and down quarks are composite and the contact interaction scale is 
^nd < 1-6 TeV or A~ d < 1.4 TeV. For a model in which all quarks are composite the 
excluded regions are A + < 1.8 TeV and A - < 1.6 TeV. 
PACS numbers: 13.87.Ce, 12.38.Qk, 12.50. Ch, 13.85.Ni 



Hard collisions between protons and antiprotons predominantly produce events 
containing two high energy jets (dijets). Measurement of the distribution of the scat- 
tering angle, between the dijet and the proton beam in the dijet center of mass frame, 
can provide a fundamental test of QCD and a sensitive probe of new physics. Di- 
jet angular distributions reflect the dynamics of the hard scattering of quarks and 
gluons, and are expected to be fairly insensitive to the momentum distributions of 
these partons within the proton. As in Rutherford scattering, dijet angular distri- 
butions from QCD processes are peaked in the forward direction. In contrast, many 
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sources of new physics produce more isotropic dijet angular distributions. Our pre- 
vious measurements using 4.2 pb _1 of data found the dijet angular distributions to 
be in good agreement with QCD predictions |IJ, and excluded a compositeness scale 
A^ d < 1.0 TeV for a contact interaction associated with compositeness of up and 
down quarks || . Here we report a measurement with a data sample that is 25 times 
larger. 

We are also motivated by the observation that the inclusive differential jet cross 
section is above a QCD prediction at high transverse energy, E T 0. Interpretations 
of this high-i?r jet excess vary from explanations within the Standard Model (mod- 
ifications of the parton distributions [|], [J or QCD corrections ||) to explanations 
beyond the Standard Model (e.g. quark compositeness 0, excited quarks 0, new Z 



bosons H , new massive gluons ||, light gluinos [ID|, and anomalous chromomagnetic 
moments of quarks [JTTj] ) . Measurements of the dijet angular distributions can help in 
resolving whether the measured excess of events with high Et jets is a signal of new 
physics or merely new information on the ingredients of QCD calculations. 

A detailed description of the Collider Detector at Fermilab (CDF) can be found 



elsewhere |12|]. We use a coordinate system with z along the proton beam direction, 
transverse coordinate perpendicular to the beam, azimuthal angle 0, polar angle 9, 
and pseudorapidity r\ = — lntan(6 l /2). Jets are identified as localized energy depo- 
sitions in the CDF calorimeters, which are constructed in a tower geometry. The 
jet axis is defined as the centroid in (77, </>) space of the calorimeter tower transverse 
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energies inside a radius R = y(A7/) 2 + (A0) 2 = 0.7 of the axis. The jet energy, E, 

— * 

and momentum, P, are defined as the scalar and vector sums, respectively, of the 
tower energies inside this radius. E and P are corrected for non-linearities in the 
calorimeter response, energy lost in uninstrumented regions and outside the cluster- 
ing cone, and energy gained from the underlying event and multiple interactions. The 
jet energy corrections increase the reconstructed jet energies on average by 24%(19%) 
for 50 GeV (500 GeV) jets. Details of jet reconstruction and jet energy corrections 



can be found elsewhere [13 . 



The dijet system consists of the two jets with the highest transverse momentum in 
the event (leading jets). We measure inclusive dijet events, defined as pp — > 2 leading 
jets + X, where X can be anything, including additional jets. The dijet invariant 



mass is defined as M = y (E± + E 2 ) 2 — (Pi + i^) 2 - We use the dijet angular variable 
X = exp(|?7i — 772D) where 771 and r] 2 are the pseudorapidities of the two leading jets. 
The variable \ has the benefit of only containing angular quantities, and hence is more 
accurately measured than a variable that involves the absolute jet energy. For the case 
of 2 — > 2 parton scattering, \ is related to the scattering angle in the center of mass 
frame, 9*, by x — (1 + 1 cos#*|)/ (1 — | cos#*|). The x distribution of QCD produced jets 
is roughly flat while many models of quark compositeness give angular distributions 
that are strongly peaked at low x- To select events with high trigger efficiency and to 
avoid problematic regions of the detector, this analysis requires x < 5, |^i| < 2, and 
1 772 1 < 2. To characterize the shape of the angular distribution in a mass bin with a 
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single number, we use the variable R x = N(x < 2.5)/iV(2.5 < x < 5), the ratio of the 
number of dijet events with \ < 2.5 to the number of dijet events with 2.5 < \ < 5. 
Isotropic angular distributions and contact interactions both tend to produce more 
events in the region \ < 2.5 than QCD, and hence will have a higher value of R x . The 
pivot point x — 2.5 was chosen to optimize the sensitivity to a left-handed contact 
interaction 0. 

Our data sample was obtained in the 1992-95 running periods using four single-jet 
triggers with thresholds on the uncorrected cluster transverse energies of 20, 50, 70, 
and 100 GeV. After applying the jet energy corrections the last three trigger samples 
were used to measure the dijet angular distribution in mass bins above 241, 300, 
and 400 GeV/c 2 , respectively. At these mass thresholds the trigger efficiencies were 
greater than 95% at all values of x considered; the average efficiency was greater 
than 99% in each mass bin. The 20 GeV trigger sample was only used to measure 
the trigger efficiency of the 50 GeV sample. The four data samples corresponded to 
integrated luminosities of 0.126, 2.84, 14.1, and 106 pb _1 . To utilize the projective 
nature of the calorimeter towers, the z position of the event vertex was required to 
be within 60 cm of the center of the detector; this cut removed 7% of the events. 
Backgrounds from cosmic rays, beam halo, and detector noise were removed with the 
cuts reported previously fL4"l , and residual backgrounds were removed by requiring 



that the total observed energy be less than 2 TeV. 

The raw x distribution was measured in five bins of dijet mass: 241 < M < 300, 
300 < M < 400, 400 < M < 517, 517 < M < 625, and M > 625 GeV/c 2 . Variations 



in the jet response and energy resolution of the calorimeter as a function of detector 77 
produced distortions in the measured angular distribution. To understand and correct 
for this effect a parametrized Monte Carlo program was developed that modeled in 
detail the measured jet response of the CDF detector after the application of the 
standard jet corrections. The relative jet response was determined from conservation 
of transverse momentum, P?, by requiring a jet in the region 0.15 < \r)\ < 0.9 and 
measuring the relative response of a jet in another pseudorapidity region (jet P T 
balancing). For jet Pt balancing, events were selected by requiring there be two and 
only two jets with Pt > 15 GeV/c, and that the azimuthal angle separating the two 
jets satisfy 150° < A<p < 210°. The largest effect was for M > 625 GeV/c 2 , where a 
6% larger jet response at \rj\ < 0.15 and a 4% smaller jet response at 0.9 < \rj\ < 1.4 
produced a tilt in the Monte Carlo angular distribution that increased the relative 
rate at x — 1 by about 10% and lowered the relative rate at x — 5 by about 10%. We 
corrected both the x an d R x distributions for these and similar effects. The correction 
reduced R x by 1%, 2%, 3%, 5%, and 6% for the 5 mass bins, respectively. The 
correction increases with dijet invariant mass because the mass spectrum is steeper 
at higher mass values, leading to a larger distortion of the angular distribution. 

In Figs. |l] and |^ and Tables I and II we present the corrected x an d R x dis- 
tributions. The data are compared to the parton level predictions of leading order 
(LO) QCD, next to leading order (NLO) QCD from the JETRAD Monte Carlo pro- 
gram pi, and QCD plus a contact interaction. For the contact interaction curve, 



we normalized LO QCD plus a contact interaction to equal the NLO QCD predic- 



tion with renormalization scale \i = Pr when the contact scale is A = oo. This was 
done by multiplying the prediction, from LO QCD plus a contact interaction, by 
the ratio of the NLO to LO QCD predictions. The LO calculations use CTEQ2L 
parton distributions, and the NLO QCD calculation uses CTEQ2M parton distribu- 
tions [16]. Alternate parton distribution sets were tried, including one in which the 



gluon distribution of the proton was significantly increased [|J, and the calculations 
were insensitive to the choice of the parton distribution. In Fig. the QCD calcula- 
tions are shown for two different choices of renormalization scale, \x = M and /i = Pt- 
The vertical axis in Fig. |^, R x , describes the shape of the angular distribution at a 
fixed mass, and is sensitive to the renormalization scale choice. The choice \x = M 
makes \i constant as a function of x i n a bin of fixed mass, while the choice // = Pp 
requires /i to vary with \. The renormalization scale dependence of the NLO calcula- 
tion, which is significantly less than that of the LO calculation, provides an estimate 
of the uncertainty in the NLO QCD calculation. Figure § also illustrates that a con- 
tact interaction would cause R x to increase at high mass, while QCD calculations 
predict that R x is roughly 0.7 at all masses shown. The angular distributions and 
angular ratio are in good agreement with the NLO QCD prediction. 

The systematic uncertainties, shown only in Fig. |2] and Table II, arise from the un- 
certainty in the jet energy response as a function of rj. The response uncertainties are 
largest in the region \rj\ < 0.15 (between 3% and 6%) and the region 0.9 < |r/| < 1.4 
(4%). Other systematic uncertainties are negligible in comparison. Since the system- 
atic uncertainties are larger than the statistical uncertainties, a cross check was per- 
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formed to verify the integrity of the measurement. The uncorrected x an d R x distribu- 
tions were remeasured with the detector pseudorapidity requirement 0.1 < \tj\ < 1.0, 
eliminating the most problematic regions of the detector. The resulting uncorrected 
distributions were then corrected back to the standard region \rj\ < 2, using the same 
parametrized Monte Carlo program, and compared with the standard results. The 
corrected x an d R x distributions for the region 0.1 < \rj\ < 1.0 agreed with those 
determined in the standard region \rj\ < 2 within the statistical uncertainties. 

The systematic uncertainties on R x are highly correlated as a function of mass. 
The diagonal terms of the covariance matrix for the R x vs. mass distribution can be 
written as Vu = af(stat)+af(sys), and the off-diagonal terms are = <7i(sys)aj(sys). 
i 7^ j, for mass bins % and j. Using this prescription and Table II the reader can recon- 
struct the full covariance matrix. We form a statistical comparison between the data 
and the theory by using the inverse of the covariance matrix, (V~ 1 )ij, and the differ- 
ence between the data and the theory in each bin, Aj, to define x 2 — Aj (V^ 1 ) a A 7 - . 
The resulting comparison between data and NLO QCD with renormalization scale 
\x = Pt is x 2 = 8-36 for 5 degrees of freedom. This is a better agreement than 
X 2 = 13.1 for NLO QCD with /i = M. 

We exclude at 95% CL any theoretical prediction which gives a x 2 of greater than 
11.1 when compared to our data. In a model of contact interactions where the up and 
down type quarks are composite we exclude at 95% CL the scales A+ d < 1.6 TeV and 
A~ d < 1.4 TeV. For flavor symmetric contact interactions among all quark flavors [117 



not just up and down quarks, the scales excluded by the angular distribution are 
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A+ < 1.8 TeV and A" < 1.6 TeV. 

We compare these exclusions with the inclusive jet cross section analysis jj, where 
we reported a broad minimum in the \ 2 between data and the compositeness model 
for the scales 1.5 < A^ d < 1.8 TeV, and best agreement for the scale A+ d = 1.6 
TeV. Here we have excluded at 95% CL the portion of this broad minimum up to the 
scale A^ d = 1.6 TeV. The inclusive jet cross section x 2 is sensitive to the choice of 
parton distributions while the dijet angular distribution \ 2 is n °t- The inclusive jet 



analysis used MRSDO' parton distributions |18j with renormalization scale \i = E T /2. 
Changing the sign of the contact interaction from positive to negative produces a 
larger deviation of the composite model from QCD for the inclusive jet cross section, 
but produces a smaller deviation from QCD in the shape of the angular distribution. 
Therefore, if we had chosen to fit the inclusive jet cross section with the negative sign 
contact interaction, best agreement would have been found with less compositeness, 
corresponding to the larger scale A~ d = 1.8 TeV, which is not excluded by the angular 
distribution. 

In conclusion, we have measured the dijet angular distributions and found them 
to be in good agreement with NLO QCD. We have presented limits on the left- 
handed contact interactions among quarks that could result if quarks were composite 
particles. Although the origin of the excess in the inclusive jet Et spectrum has not 
been determined, the angular distribution data are consistent with the hypothesis 
that the high-.&r jet excess is caused by effects within the Standard Model. The 
angular distribution data exclude at the 95% CL the hypothesis that the high-.&r 
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jet excess is caused by a contact interaction among up and down quarks with scale 
A+ < 1.6 TeV. 
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Table I: The dijet angular distribution and statistical uncertainty for the five mass 
bins (GeV/c 2 ). 







(100/N)(dN/dx) 






X 


241 < M < 300 


300 < M < 400 


400 < M < 517 


517 < M < 625 


M > 625 


1.25 


31.1 ±0.7 


31.7±0.5 


31.9 ±0.5 


32.6 ± 1.2 


31.7 ±2.4 


1.75 


26.8 ±0.6 


26.5 ±0.5 


26.3 ±0.4 


27.2 ± 1.1 


26.5 ±2.2 


2.25 


23.0 ±0.6 


23.8 ±0.5 


24.3 ±0.4 


25.1 ± 1.1 


26.3 ±2.2 


2.75 


23.4 ±0.6 


23.2 ±0.5 


23.9 ±0.4 


23.0 ± 1.0 


25.3 ±2.2 


3.25 


24.3 ±0.7 


23.8 ±0.5 


23.5 ±0.4 


21.2 ± 1.0 


21.6 ±2.0 


3.75 


22.5 ±0.6 


24.0 ±0.5 


23.3 ±0.4 


24.1 ± 1.1 


22.2 ±2.1 


4.25 


24.6 ±0.8 


23.3 ±0.5 


23.6 ±0.5 


22.8 ± 1.1 


22.4 ±2.1 


4.75 


24.4 ±0.8 


23.7 ±0.5 


23.1 ±0.5 


24.1 ± 1.2 


24.0 ±2.2 



Table II: The mean dijet mass (GeV/c 2 ), number of events, dijet angular ratio R x , 
and its statistical and systematic uncertainty. The completely correlated systematic 
uncertainty can be used to form the covariance matrix (see text). 



<Mass> 


Events 


R x 


Stat. 


Sys. 


263 


15023 


0.678 


0.012 


0.018 


334 


23227 


0.695 


0.010 


0.025 


440 


28202 


0.703 


0.009 


0.033 


557 


4425 


0.738 


0.023 


0.054 


698 


1056 


0.732 


0.046 


0.103 
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Figure 1: The dijet angular distribution (points) compared to predictions of NLO 
QCD (solid curve), LO QCD (dashed curve), and LO QCD with a quark contact 
interaction (dotted curve). The contact interaction calculation is normalized to equal 
NLO QCD when A u d = oo (see text). Error bars on the data and NLO QCD are 
statistical. 
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Figure 2: The dijet angular ratio (points) as a function of the dijet invariant mass, 
compared to LO QCD (dashed curve), NLO QCD (solid curve), and LO QCD with 
a quark contact interaction normalized to NLO at A ud = oo (dotted curve). QCD 
is shown for two renormalization scales (/i = M and \x = Pt)- Contact interactions 
are displayed for three different compositeness scales, with two different signs for the 
amplitude of the contact term (upper dotted curve is A+ d , lower dotted curve is A~ d ). 
The inner error bars on the data are statistical uncertainties and the outer error bars 
are statistical and systematic uncertainties added in quadrature. The error bars on 
NLO QCD are statistical. 
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